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We derive general bounds on the linear response energy absorption rates of periodically driven 
many-body systems of spins or fermions on a lattice. We show that for systems with local interac¬ 
tions, energy absorption rate decays exponentially as a function of driving frequency in any number 
of spatial dimensions. These results imply that topological many-body states in periodically driven 
systems, although generally metastable, can have very long lifetimes. We discuss applications to 
other problems, including decay of highly energetic excitations in cold atomic and solid-state sys¬ 
tems. 
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INTRODUCTION 

Time-dependent driving recently emerged as a new 
versatile tool for engineering various quantum states of 
matter. In non-interacting systems, periodic driving can 
be used to modify band structures, and in particular to 
make them topologically non-trivial[l-4]. 

Experimentally, periodic driving has been used to re¬ 
alize strong artificial magnetic fields [5], as well as 2D 
Bloch bands with non-zero Chern numbers (similar to 
the Haldane model) in systems of cold atoms in opti¬ 
cal lattices [6, 7]. Theoretically, natural extensions of 
these ideas to topologically non-trivial Floquet many- 
body states have been proposed. In particular, it was 
suggested that topological Floquet bands can host frac¬ 
tional Chern insulators [8] , as well as symmetry-protected 
topological states [9]. 

Topological order and more generally quantum order 
are usually associated with ground states and low tem¬ 
peratures, where they are protected by a finite excitation 
gap. However, periodic driving breaks energy conserva¬ 
tion, making the very concept of the ground state mean¬ 
ingless. It was argued [16-18] that under driving, generic 
many-body systems that obey the eigenstate thermal- 
ization hypothesis [10-12], eventually heat to up to an 
infinite-temperature, featureless state. Thus, ’’Floquet 
topological insulators” are generally metastable. It is 
important to understand their lifetimes, and use the the¬ 
oretical understanding to design experiments in which 
the Floquet many-body states would be long-lived. 

In this paper, we derive general results regarding the 
heating of periodically driven many-body systems on a 
lattice. We consider both the cases of local driving (time- 
dependent perturbation acting only on a few degrees of 
freedom), and a global driving (driving applied every¬ 
where in the system). The latter setup is relevant to cold 
atoms experiments [5-7]. Assuming that interactions are 
local, we prove a general bound for the linear-response 


heating rates, which indicates that at high driving fre¬ 
quency (much higher than a natural energy scale of the 
system, e.g., kinetic or interaction energy of one parti¬ 
cle), heating is exponentially slow. Fundamentally, this 
bound follows from the locality of quantum dynamics in 
systems with local interactions, and, for the case of global 
driving, relies on the Lieb-Robinson bounds. 

RESULTS 

We consider a lattice system of spins or fermions 
with a local Hamiltonian H = subject to a 

periodic time-dependent perturbation with an operator 
O = Oi, which is a sum of one or more local terms Oi. 
It is understood that i runs over the sites of the lattice 
and hi, Oi act on a fixed, finite number of sites around i. 
For simplicity, we focus on the case of harmonic driving 
with frequency ui and strength g: 

H{t) = H-\-gcos{(jjt)0. (1) 

and we fix 11 Oil I < 1 for concreteness. We assume 
that the system is initially in a thermal equilibrium, 
and is described by a density matrix pp, /? = 1/T. 
The energy absorption rate dE/dt is related to 
the dissipative part of the linear-response function, by 
dE/dt = 2g'^uia{uj), upon averaging over a cycle and up 
to higher orders in g. The response cr(a;) can be expressed 
in terms of different-time commutators of operators Op. 

,, ^ d-oo 

( 2 ) 

where (...)/? denotes thermal averaging and Oi{t) refers 
to the Heisenberg dynamics generated by the time- 
independent H. One can show that uja(uj) is positive 
and symmetric in w, see e.g. [32, 33]. For concreteness, 
let us take w > 0 so that cr(uj) > 0. To be precise, in fi¬ 
nite volume, uja(oj) is a positive distribution rather than 
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a bonafide function and we will need to integrate this 
distribution over small intervals to state rigorous results: 
we write /([wi, W 2 ]) = /(w) duj. 


Local driving 

First, using the locality of the Hamiltonian, we prove 
the following bound for individual terms aij(uj): 

There is a k > 0 and a numerical constant C such 
that, for any 6oj > 0,oj > 0, 

|ay([w,cc + H)l (3) 

uniformly in the total volume. Note that it is only the dis¬ 
sipative (asymmetric in time) part of the response func¬ 
tion that decays exponentially. The reactive (symmetric 
in time) part contains a principal value integral rather 
than a J-function, cfr. (5), in energy, and it generically 
decays no faster than l/w. 

We first prove (3) for i = j. It is convenient to write 
an as a sum of contribution of individual eigenstates |n) 
oi H: 


au(uj) = 7r^p„(7”(a;) - 7D(-w)), (4) 

n 


where ^ ^ , Z = Tr(e is the probability that 

the system is in eigenstate |n), and 7 )) is given by 

= 51 \{'m\0,\n)\^6{Er, - - uj). (5) 


To estimate this quantity, let us rewrite it as follows, for 
any k: 


» =i: 




CJ 


2k 


d{En Ejji io), 

( 6 ) 


where the r.-h.s. contains k commutators with iJ, and we 
have used the fact that En — Em = oj. Next, we use the 
fact that in local systems, the norm of [[[O^, H],E[],El] 
can be bounded as 


\\[[[0,,H],H],...,H]\\<e’^k\, (7) 

where e is an energy scale that can be expressed via 11 11 

(local norm of the Hamiltonian), the range of hj and Oi, 
and the coordination number of the lattice. 

Integrating Eq.( 6 ) over an interval [u},u} + 6uj\, and us¬ 
ing inequality (7), we obtain: 

7))([u;,a; + H)< (^) < (^) -VfceN (8) 

Choosing k = jj, we arrive at the estimate 

7))([w,a; + H)<e-"“, «=-, (9) 

ee 


and therefore, using (4) and the fact that J2nPn = 1, we 
obtain 


aii{[uj,u + 5u\) <2TTe (10) 


The off-diagonal terms Uij, i j, can be bounded by 
diagonal ones using positive-definiteness of ujaij and the 
Cauchy-Schwartz inequality: 


\a^j{uj)\ < 


1 


{au{u}) + ajj{u})). 


( 11 ) 


Thus, just like diagonal terms, off-diagonal terms decay 
exponentially at large frequency, concluding the proof. 

The above result can be immediately applied to the 
case of local driving, when the operator O acts only on 
a finite number of n lattice sites (i.e. Oi on other sites 
are taken to be zero), while the system size N is taken 
to infinity. For such a setup, the energy absorption rate 
will be smaller than ~ where n is the number 

of lattice sites affected by the periodic driving. 

The result (3) has a clear intuitive meaning for two par¬ 
ticular cases. The first example is a gapped many-body 
system at zero temperature, e.g., a gapped spin chain 
with a maximum energy per spin of order J. The exis¬ 
tence of the gap implies exponential decay of correlations 
with a typical length scale When the spin chain is sub¬ 
ject to a periodic drive, in order to absorb a large energy 
w, a local operator Oi must flip a large number n ^ uj/ J 
spins. It is clear that at large n ^ f the amplitude of 
such a process should be exponentially small. Another 
example is provided by a system with weak interactions, 
e.g., a weakly interacting Fermi-gas. Then, absorbing a 
large energy w requires a creation of n ~ tu/J particle- 
hole pairs (J in this case denotes a maximum energy of 
one electron-hole pair) - a process which will be sup¬ 
pressed as H", where V is the interaction strength. Our 
analysis shows that the energy absorption rate remains 
exponentially small much more generally: the bound ap¬ 
plies to systems without gaps and also to systems with 
strong interactions. 


Global driving 

Next, we consider the case of global driving, when Oi 7 ^ 
0 for all lattice sites. This setup is generic and relevant, 
in particular, to all experiments in which driving is used 
to create topological states. In this case, we prove a 
different bound: 

There is a k > 0 and a numerical constant C such 
that, for any 6 lj > 0,uj > 0, 

|(T([a;,w-b(jw])| < 7VC'e-'=l“l. (12) 

with N the total number of spins. It is intuitively clear 
that the absorption rate in this case should be propor¬ 
tional to N. Therefore one cannot simply use the bound 
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(3), since it would give a heating rate that scales as 
N^. Thus, we use Lieb-Robinson bounds [38] to esti¬ 
mate Uij for remote i and j. For two operators A,B 
with ||A||,||i?|| < 1 and support in regions X,Y, they 
read: 

\\[A{t),B]\\<Ce-<^-^--^\ r = dist{X,Y), (13) 

with vlr the Lieb-Robinson velocity and C a numerical 
constant (below we always use C for numerical constants 
whose value can change from line to line). We choose 
to measure distances in units of the lattice spacing, and 
hence a is simply a numerical constant as well, a, C, vlr 
are determined by the norm and range of hi and the 
type of lattice. Without loss of generality, we assume 
KUJ ^ 1 (for small values of kw, the bound can always be 
satisfied by tuning C). We will prove the bound (12) for 
5ui = (5a;)o with some arbitrary but fixed {Suj)o > 0 (we 



uj'—uj \ ^ 

) o - ij ( a ;') = ^/tt{6uj)o / 

^ —OO 


Our strategy is to estimate terms with dist(i,j) > r*, 
with r* large, using Lieb-Robinson bounds, and to bound 
terms with dist(i,j) < r* using (3). We will choose r* = 
2kujla with k being the k, featuring in (3). 

Fist, we study a general (i,j)-term in (15) with 
dist(*,j) = r > r^. We break the time-integral in the 
r.h.s. of (15) into an integral over the interval [—tc,tc] 
with tc = tc(r) = r/(2vLii), and an integral over the rest 
of the real axis. The former is bounded by LR-bounds 
(13): 




'-tc 




(16) 

and the latter is bounded using ||[Oi(t),Oj]|| < 2 (since 

lia||<l): 
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(17) 

Using (16,17) and setting St < tc (since {S(jj)o is fixed and 
KUJ ^ 1 ), we bound the sum over (z, j) with r > r* by 


CNr^-^ 


-arj2 

Vlr 



(18) 


By increasing r*, the exponent in the second term be¬ 
comes at least as small as that in the first term and we 
can bound (18) by 


CNrf 





will use in the proof that (<5a;)o is smaller than quantities 
that diverge with ui —> oo). The bound for arbitrary Su 
is then recovered as follows: For Suj < ((5a;)o, we have 
CT([a;, w-I- (5cli]) < ct([w, w-|- ((5w)o]), hence the bound. For 
Suj > (Suj)o, we dominate a{[uj,uj + ^wj) < 
fc((5a;)o,a;-|-(fc-|-l)(i5a;)o]) and we apply the result for each 
term in the sum, so that the result follows by readjusting 
C. 

First, we dominate, for uj > 2(Suj)o, 


cr([a;,w-k ((5a;)o]) < - 

1 — e 



e 



a(uj') 


(14) 

This relies on positivity of CT(a; > 0) and the symmetry 
a{—uj) = —ct(uj). Then we split cr(uj) — in 

(14) and we recast the resulting integrals in the time 
domain: 




St = 


(Suj)o' 


(15) 


This provides a bound for the contribution of remote 
pairs \i — j\ > r* to the r.-h.s. of Eq.(15). Multiplying by 
y/TT{Suj)(i and using {Suj)q/vlr < C, we bound their con¬ 
tribution to the response function (14) by CNr^e~'^‘^. 
There are ^ rf remaining terms with dist(z,j) < r*. 
Their contribution to (14) can be bounded using (3), 
which gives CNr^e~'^‘^, as well. Recalling r* = 2kujja, 
we get an overall bound (i.e. summed over all r) of the 
form 


CN{kujfe-~^^ 

for the response function. By slightly reducing k (the 
new value is called k again) and increasing C, we get the 
bound ( 12 ). 

The necessity of averaging over a frequency window Suj 
is likely an artefact of our proof. We expect that in the 
thermodynamic limit, the bounds (3,12) hold for a(uj) 
itself, i.e. that the distribution is a bonafide function, but 
we cannot prove this, see also [36, 37] for mathematical 
details and [35] for a polynomial bound at large uj. 


EXTENSIONS 

Let us briefly comment on extensions of our results. 
First, we note that in d = 1 we can in fact choose the de¬ 
cay rate k to he arbitrarily large at the cost of increasing 
the numerical prefactors C in (3,12). This is achieved by 
improving the bound in the r.h.s. of (7) to C{'y)e~~^'^kl, 
for any 7 > 0, as described in [28, 39]. The extension 
for local driving is immediate and for global driving, we 
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then simply choose r* with a larger k. Hence in d = 1, 
the decay is in fact superexponential, but one should not 
expect such an improvement to hold in d > 2 [29], nor 
for the case of quasilocal, exponentially decaying Oi^hi 
(instead of strictly local). 

It is also possible to obtain similar results as ours for 
some models of lattice bosons, e.g. the Bose-Hubbard 
model at high temperature, but then the exponential de¬ 
cay is weakened to a stretched exponential. The key ob¬ 
servations are that a) the bound (7) fails trivially for un¬ 
bounded hi and one has to use a weighted norm instead 
and h) one needs to assume (spatial) decay of correlations 
for the thermal ensemble {■)p (provable by cluster expan¬ 
sions at high temperature) as Lieb-Robinson bounds are 
no longer available, see [27]. 

The most important extension [27], however, is to go 
beyond linear response and to allow for general initial 
states, showing that the phenomenon described here is 
quite similar to localization in energy^ except that, pre¬ 
sumably, it in general breaks down after a sufficiently 
long time, [16-18]. 

Finally, for completeness, we mention that there are 
two remarkable cases where the energy localization does 
not break down for long times: 1) driven MBL systems 
at not too small frequencies, see [21-23] and 2) non¬ 
interacting fermions, i.e. with hi,Oi containing only lin¬ 
ear and quadratic terms in c, cb In that case, by per¬ 
forming a canonical transformation, we fall back on a 
one-particle problem where exact dynamical localization 
is possible [34]. In the case where H describes non¬ 
interacting fermions but O is of order q > 2 in the fields 
c, , we do not expect genuine localization but still the 
linear response vanishes exactly: cr(w) = 0 whenever w 
exceeds the bandwidth times q, as one sees from (5). 

DISCUSSION 

The main implication of the above results is that, 
although many-body states (e.g. topological states) 
in isolated, periodically driven systems are generally 
metastable, they have a very long life time, if driving 
frequency is much higher than the natural energy scale 
of the system. We note that this limit is indeed realized 
in recent experiments. 

Further, we note that there has been recent interest in 
dynamical localization in periodically driven many-body 
systems [17, 19-24]. In such studies, numerical simula¬ 
tions are a useful tool. Our results imply that, in order 
to observe delocalization at high driving frequency, one 
may have to study the system dynamics at (exponen¬ 
tially) long times. 

Our results also can be directly applied to a different 
class of problems: the decay of a highly energetic exci¬ 
tation into many low-energy excitations. One physical 
model where such a problem naturally arises and has 


been studied is the large-!/ Fermi-Hubbard model, 

H = J ^ cl^Cjs + U'^ni^nii, U :$> J (19) 
(b>,s=t,4. » 

where (ij) denotes neighbouring sites, and s is a spin 
label. Doublons (doubly occupied sites) have a typical 
energy ~ U which is much greater than the kinetic energy 
J. At temperatures T U, when there are very few 
doublons in the system, one can ask how quickly doublons 
decay into particle-hole pairs. This problem has been 
addressed experimentally [25] and theoretically [26] using 
perturbation theory, and the decay rates were found to 
be exponentially small in U/J. Therefore, in this case 
our bound appears to be saturated. 

We note that the bound is non-perturbative in the in¬ 
teraction strength of the systems, and we expect it to be 
useful for other strongly interacting systems where exci¬ 
tations with very different energy scales are present, e.g., 
equilibration of Fermi-Fermi cold gases with very differ¬ 
ent mass parameters, as well as random spin models with 
a broad distribution of exchange couplings. For random 
spin models, an exponential decay of a local spin cor¬ 
relation function at high frequency has been previously 
obtained using mean-field-type approximations [30, 31]. 
Our results are consistent with this earlier study in d > 1, 
but show that in one-dimensional systems the decay at 
high frequency is even faster than exponential. 

CONCLUSION 

We have proven that the dissipative part of the linear 
response for local lattice systems decays exponentially 
at high frequencies. In particular, this means that 
heating by periodic driving will be exponentially slow. 
This result provides a foundation for so-called Floquet 
many-body phases by showing that, though metastable, 
they will be very long-lived. 
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